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Abstract 

A confluence of numerical and theoretical results leads us to conjecture that the Hilbert-Schmidt 
separability probabilities of the 15- and 9-dimensional convex sets of complex and real two-qubit 
states (representable by 4 x 4 density matrices p) are ^ and respectively. Central to our 
reasoning are the modifications of two ansatze, recently advanced {Phys. Rev. A, 75 [2007], 
032326), involving incomplete beta functions B v (a,b), where v = p^lptt ' now ' se ^ the 

separability function S rea i(u) proportional to B v {y, ^,2) = |(3 — Then, in the complex 

case — conforming to a pattern we find, manifesting the Dyson indices (/3 = 1,2,4) of random 
matrix theory-we take S comp i ex {v) proportional to S^ eal (v). We also investigate the real and 
complex qubit-qutrit cases. Now, there are two variables, v\ = , v<i = , but they 

appear to remarkably coalesce into the product n = v\Vi = , so that the real and complex 

separability functions are again univariate in nature. 
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I. INTRODUCTION 

Zyczkowski and Sommers have derived — using random matrix theory (in particular, the 
Laguerre ensemble) — general formulas for the (n 2 — l)-dimensional and the 
n(n+i)-2 _ c [j mens j ona i volumes of the complex and real nxn density matrices (p), respectively, 
in terms of the Hilbert-Schmidt (HS) metric , sec. 14.3] (as well as the Bures metric 
0, sec. 14.4] 3, J]). Later, Andai jq] examined these and related questions, using a quite 
different framework. He applied mathematical induction on the leading principal minors 
of p, along with the established formulas for hyperareas of surfaces of n-spheres and beta 
integrals. He reproduced — up to normalization factors — the HS real and complex volume 
formulas in [1] (and, moreover, the (2n 2 — n — l)-dimensional quaternionic volumes). (In 
addition to the HS and Bures metrics, Andai considered, for the single qubit case, the broad 
[infinitely nondenumerable] class — which does include the Bures as its minimal member — of 
monotone metrics. Unlike Zyczkowski and Sommers, he did not obtain formulas for the 
hyperareas occupied by density matrices of less than full rank.) 

Despite these considerable theoretical advances, volume (and, hence, probability) for- 
mulas have not yet become available for the important subsets of separable (n < 6) and 
positive-partial-transpose (n > 8) nxn density matrices (n composite). (Szarek [6( employed 
methods of asymptotic convex geometry to estimate the volume of the set of separable mixed 
quantum states for N qubits, and Aubrun and Szarek [?J for iV qudits. It was concluded 
in these studies that the separable volumes were superexponentially small in the dimension 
of the set of states. For large D, the (D 2 — l)-dimensional volume for bipartite systems 
of positive-partial-transpose states, however, is much larger than the volume of separable 
states 17, Thm. 41.) 

n 

To address this fundamental lacuna, at least in the Hilbert-Schmidt context (cf. [8|), 



we developed in [9] a methodology — incorporating the Bloore parameterization of density 
matrices [lOj (sec. II Ah . Its numerical application led to ansatze, involving (apparently inde- 
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pendent) incomplete beta functions, for the 9- dimensional real and 15-dimensional complex 
separable volumes in the qubit-qubit (n = 4) case [9|. In the sequel to that study here, 
we, first, apply this Bloore framework to various scenarios involving n x n density matri- 
ces (n = 4,6,8,9), in which certain of their off-diagonal entries have been nullified. This 
enables us to now obtain exact results, of interest in themselves, and possibly suggestive 
of solutions/approaches to the full (non-nullified) highly computationally-challenging prob- 
lems. 

In fact, based on certain (real-complex-quaternionic) patterns emerging in these exact 
results (sec. E 



matrix theory 



), bearing an obvious relation to the Dyson indices ((3 = 1, 2, 4) of random 
111 ], we are led to modify the incomplete beta function ansatze for the two full 
(real and complex) problems advanced in [9j. The "separability function" in the complex 
case is now not analyzed as if it were independent of that in the real case (which we still take 
to be an incomplete — but slightly different — beta function), but actually simply proportional 
to its square. These central analyses will be elaborated upon in sec. HXl (eqs. (194]) - ( l96l) ). 
where it shown that the modified ansatze do, in fact, accord well (Fig. [3]) with the numerical 

LI 

results of [9J. 

We begin our extensive series of lower-dimensional analyses, by examining a number of 
two-qubit scenarios (sec. [IT])- In them, we are able to compute a number of interesting 
exact two-qubit scenario-specific HS separability probabilities. (Listing them in increasing 
order, we have {i §, §, f , ±§f , ^, ff , §, ^gf , 2 - f§f , ±|, l}.) For each of the scenarios, we 
identify a certain univariate separability function S scenar i (is), where v = . The integral 
over v G [0, oo] of the product of this function (typically of a piecewise nature over [0,1] 
and [l,oo]) with a scenario-specific (marginal) jacobian function J sce nario{y) yields the HS 
separable volume iV^~). The ratio of V^f to the HS total (entangled and non-entangled) 
volume (V t o t s ) gives us the HS scenario-specific separability probability. 

The question of the "relative proportion" of entangled and non-entangled states in a given 
generic class of composite quantum systems, had apparently first been raised by Zyczkowski, 
Horodecki, Sanpera and Lewenstein (ZHSL) in a much-cited paper jl2| . They gave "three 
main reasons" — "philosophical", "practical" and "physical" — upon which they expanded, 
for pursuing the topic. The present author, motivated by the ZHSL paper, has investigated 
this issue in a number of settings, using various (monotone and non-monotone) measures on 
quantum states, and a variety of numerical and analytical methods {4, l], 14, 3, 16, 17, 3, 
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19j (cf. [20 



21 



22 



23j). Though the problems are challenging (high-dimensional) in nature, 



many of the results obtained in answer to the ZHSL question in these various contexts have 
been strikingly simple and elegant (and/or conjecturally so). 

Specifically here, we further develop the (Bloore-parameterization-based) approach pre- 
sented in [9]. This was found to be relatively effective in studying the question posed by 
ZHSL, in the context of two-qubit systems (the smallest possible example exhibiting entan- 
glement), endowed with the (non- monotone 24]) Hilbert- Schmidt (HS) measure [lL inducing 
the flat, Euclidean geometry on the space of 4x4 density matrices. This approach [9J exploits 
wo distinct features of a form of density matrix parameterization first discussed by Bloore 
1(3 ]. These properties allow us to deal with lower- dimensional integrations (more amenable 
to computation) than would otherwise be possible. We further find that the interesting 
advantages of the Bloore parameterization do, in fact, carry over — in a somewhat modified 
fashion — to the qubit-qutrit (sec. Illlland IX Al) . qutrit-qutrit (sec. HVj) and qubit-qubit-qubit 
(sees. |V] and IVT]) domains. 



A. Bloore (off-diagonal-scaling) parameterization 



We, first, consider the 9-dimensional convex set of (two-qubit) 4x4 density matrices with 
real entries, and parameterize them — following Bloore [l(| (cf. 0, P- 235]) — as 

( Pll Z\1y/P\\Pl1 Zi^pup^ Zu^PllPM \ 



P 



Z\2\/ ' P\\Pl2 P22 ^23VP22P33 ^24a/P22P44 

z 13y/PllP33 Z 23y/P22p33 P33 %4v / P33P44 
\ Zi4y/pnpu Z24y/P22Pu ^34 ^ P33P44 P44 



(1) 



One, of course, has the standard requirements that pa > and (the unit trace condition) 
YiiPa = 1. Now, three additional necessary conditions (which can be expressed with out using 
the diagonal entries, due to the pa > stipulation) that must be fulfilled for p to be a density 
matrix (with all eigenvalues non-negative) are: (1) the non-negativity of the determinant 
(the principal 4x4 minor), 



(^34 — l) ^12 + 2 (^14 (^24 — ^23^34) + ^13 (^23 — ^24^34)) ^12 



-23 ^24 



(2) 



+ ^14 { Z 23 _ l) + Z 13 { Z 24 - l) + 2^23^24^34 + ^Z 13 Z 1A (z U - Z23Z24) + 1 > 0; 
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(2): the non-negativity of the leading principal 3x3 minor, 

— z\ 2 + 2^13^23^12 — z\ 3 — + I > 0; 



(3) 



and (3): the non-negativity of the principal 2x2 minors (although actually only the i = j = 1 
case is needed, it is natural to impose them all), 

1 - 4 > 0. (4) 

As noted, the diagonal entries of p do not enter into any of these constraints — which taken 
together are sufficient to guarantee the nonnegativity of p itself — as they can be shown to 
contribute only (cancellable) non-negative factors to the determinant and principal minors. 
This cancellation property is certainly a principal virtue of the Bloore parameterization, 
allowing one to proceed analytically in lower dimensions than one might initially surmise. 
(Let us note that, utilizing this parameterization, we have been able to establish a recent 
conjecture of Mansson, Porta Mana and Bjork regarding Bayesian state assignment for three- 
evel quantum systems, and, in fact, verify our own four-level analogue of their conjecture 



2J, eq. (52)] |27j.) 



28. 



29|, 



30j requiring the non- 



Additionally, implementing the Peres-Horodecki condition 
negativity of the partial transposition of p, we have the necessary and sufficient condition 
for the separability (non-entanglement) of p that (4): 

v (zg 4 - l) z\ 2 + 2y/u (uz 13 z u + z 23 z 2i - \fv {z u z 23 + z 13 z 24 ) z u ) z l2 - z\ z - vz\± + u+ (5) 

+V {{z 2 2A - l) zf 3 - 2Z U Z 23 Z 24 Z 13 - zl 4 + Z 2 U [z\ z - v)) + 2^(^32:23 + ^14^24) ^34 > 0, 

where 

2 P11P44 la , 
v = p = , (6) 

P22P33 

being the only information needed, at this stage, concerning the diagonal entries of p. (It 
is interesting to contrast the role of our variable v, as it pertains to the determination 
of entanglement, with the rather different roles played by the concurrence and negativity 



31 



32].) We have vacillated between the use of v and p as our principal variable in our two 
previous studies 0, Q. In sec. IVIH we will revert to the use of p, as it appears that its use 
can avoid the appearances of square roots, which, it is our impression, at least, can impede 
certain Mathematica computations. 
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1. Reduction of dimensionality 



Thus, the Bloore parameterization is evidently even further convenient here, in reducing 
the apparent dimensionality of the separable volume problem. That is, we now have to 
essentially consider only the separability variable v rather than three independent (variable) 
diagonal entries. (This supplementary feature had not been commented upon by Bloore, 
as he discussed only 2x2 and 3x3 density matrices, and also, obviously, since the Peres- 
Horodecki separability condition had not yet been formulated in 1976.) The (two variable — 
Ui,v 2 ) analogue of (jSJ) in the 6x6 (qubit- qutrit) case will be discussed and implemented 
in sec. IHII Additionally still, we find a /bur-variable counterpart in the 9x9 qutrit-qutrit 
instance [sec. IIVJ . and three- variable counterparts in two sets of qubit-qubit-qubit analyses 
(sees. |V] and IVIj) . (The question of whether any or all of these several ratio variables are 
themselves observables would seem to be of some interest.) It certainly appears to us that in 
the qubit-qutrit case (sec. IHII and IX Al) the two associated ratio variables (p±, v 2 ) importantly 
merge or coalesce into the simple product r) — v\v 2 for all analytical purposes. (Products of 
ratio variable do also appear in the limited number of still higher-dimensional analyses we 
report below, so perhaps some similar merging or coalescing takes place in those settings, 
as wehV) 

eqs. (3)-(5)], we expressed the conditions (found through application of the "cylin- 



In 



drical algebraic decomposition" 34|) that — in terms of the Bloore variables arbi- 



trary 9-dimensional 4x4 real density matrix p must fulfill. These took the form, 

2 12i ^13) Z U £ [ — 1) 1]) ^23 £ [-^23' ^23]) ^24 £ [^24) z U [^34> ^34]) (?) 



where 



Zfz = z 12 z 13 ± y 1 - zf 2 yj 1 - zf 3 , Z 24 = z 12 z 14 ± y 1 - z\ 2 \j 1 - zj 4 , (8) 

7 ± _ Z13Z14 — Z\ 2 Z\ 4 Z 2 -& — Z12Z13Z2A + ^23 ^24 i S 



^34 



-12 



and 



1 + Z 12 + Z 13 — 1Z\ 2 Z\ 3 Z 23 + Z 23 \J —I + Z 12 + z i4 — 2212^14^24 + Z 24 - (9) 
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2. Possible transformations of the 's 



In his noteworthy paper, Bloore also presented [lOj, sees. 6,7] a quite interesting discussion 
of the "spheroidal" geometry induced by his parameterization. This strongly suggests that it 
might prove useful to reparameterize the Zij variables in terms of spheroidal-type coordinates. 
Following the argument of Bloore — that is, performing rotations of the (2:13, Z23) and (214, Z24) 
vectors by j and recognizing that each pair of so-transformed variables lay in ellipses with 



axes of length \/l ± z 12 — we were able to substantially simply the forms of the feasibility 
conditions (©-©). 

Using the set of transformations (having a jacobian equal to ^ ) 



y v *12 ouo \u\ j -r 0111 \y\ j v ^12 t v /l /2 ~r j- /lnN 
z iz — > 7= j C-LU; 



^23 



^14 



^24 





- Z\2 cos 


A) + sin(^) 


V z 12 + 1 j V7l72 + 1 


(sin ( 








\Ai2 - 


f 1 - COS (0i ) 


V 1 - 212 j V7l72 + 1 






y/2 






— Z\2 COS 


(0 2 ) + sin(0 2 ; 


V«12 + l) \/7l + 72 


(sin 






^2) >/*12 


+ 1 - COS (0 2 , 


Vl - 212 J V7l + 72 



cos (0i - 6 2 ) V71 + 72 V7i72 + 1 
£34 — ► ^34 ^= j 

one is able to replace the conditions (©-(El)) that the real two-qubit density matrix p — given 
by ©—must fulfill by 

71 e [0,1]; 72 e[7i,-]; Z 34 e [-71,71]; z 12 e[-l,l]; 0i,0 2 e [0,2tt]. (11) 
7i 

We developed this set of transformations at a rather late stage of the research reported here, 
and presently have no indications that are of any special aid in regard to the particular 
difficulties/challenges posed by the HS separability-probability question. So, the qubit-qubit 
results reported below (sec. [II]) do rely essentially upon the conditions (©-([9])) and the 
original parameterization in terms of the z^s of Bloore. 

Another very interesting simplifying parameterization — expressed in terms of correlations 



and partial correlations — can be found in the statistical/mathematical literature 25 



37 



35, 



36, 



. (In fact, the Bloore parameterization can be readily seen — in retrospect — to be 



simply a way of decomposing a density matrix into a correlation matrix (cf. |39|), plus its 

9 



diagonal entries.) But this too seems to have no particular enhanced value in analyzing 
partial transposes. (The cited literature also appears to be highly relevant to the problem 
of the random generation of density matrices.) 



B. Previous analysis and beta function ansatze 



Is 



In |9j, we studied the four nonnegativity conditions (as well as their counterparts — having 
completely parallel cancellation and univariate function properties — in the 15-dimensional 
case of 4 x 4 density matrices with, in general, complex entries) using numerical (primarily 
quasi- Monte Carlo integration) methods. We found a close fit to the function [9, Figs. 3, 4], 
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-•approx 1 
real ' 



( 4+ ^) B fr^) *(?^)- < i2 » 



entering into our formula (cf. ( 1TT1) . ( fi~8i) ' 

'1 POO 



K^/real = 2 / Jreal{y)S rea l[y)dv = / J real (v)S rea i(v)dv, (13) 

Jo Jo 
for the 9-dimensional Hilbert-Schmidt separable volume of the real 4x4 density matrices 9j, 
eq. (9)]. Here, B denotes the (complete) beta function, and B v the incomplete beta function 



B v (a,b) = / w a - l {l - w) b - l dw. (14) 
Jo 

Additionally |9|, eq. (10)], 

_ v*l 2 (12 (u(u + 2) (y 2 + Up + 8) + 1) log (y/u) - 5 (5^ 4 + 32z/ 3 - 32u - 5)) 

JreaiW- 3780(^-1)9 ^) 

is the (apparently highly oscillatory near v = 1 [9_|, Fig. 1]) jacobian function resulting from 
the transformation to the v variable of the Bloore jacobian (Ii-f = iPu) % ■ (A referee did indicate 
that the apparent oscillations vanished, when he employed a Maple program using 50 digits 
of precision. [Only in the latest Version 6 of Mathematica is a comparable plot feasible.] 
Also, perhaps, we should refer to J rea i[y) as a marginal jacobian, since it is the result of the 
integration of a i/iree-dimensional jacobian function over two, say pn and P22, variables.) 
In the 15-dimensional complex two-qubit case, we found that the function 

/ 100000000 \ /2V6 3 V 4 /V6 3\ 

^»=(^fTifJ b~'^J "[—'T2)' (16) 

provided a close fit to our numerical results [18|, eq. (14)]. 
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C. Research design and objectives 



Although we were able to implement the three (six-variable) nonnegativity conditions 
(©, © and (|1|)) exactly in Mathematica in for density matrices of the form (pQ), we 
found that additionally incorporating the fourth Peres-Horodecki (separability) one (151)— 
even holding v fixed at specific values — seemed to yield a computationally intractable prob- 
lem. 

In light of the apparent computational intractability in obtaining exact results in the 
9-dimensional real (and a fortiori 15-dimensional complex) two-qubit cases, we adjusted the 
research program pursued in We now sought to determine how far we would have to 
curtail the dimension (the number of free parameters) of the two-qubit systems in order to be 
able to obtain exact results using the same basic investigative framework. Such results — in 
addition to their own intrinsic interest — might help us understand those previously obtained 
(basically numerically) in the full 9-dimensional real and 15-dimensional complex cases q| 
(which, retrospectively, in fact, we assert does turn out to be the case). 

To pursue this lower-dimensional exact strategem, we nullified various m-subsets of the 
six symmetrically-located off-diagonal pairs in the 9-parameter real density matrix (III) , and 
tried to exactly implement the so- reduced non-negativity conditions ((J2D, (J3|), (jlj and ©)— 
both the first three (to obtain HS total volumes) and then all four jointly (to obtain HS 
separable volumes). We leave the four diagonal entries themselves alone in all our analyses, 
so if we nullify m pairs of symmetically-located off-diagonal entries, we are left in a (9-m)- 
dimensional setting. We consider the various combinatorially distinct scenarios individually, 
though it would appear that we also could have grouped them into classes of scenarios 
equivalent under local operations, and simply analyzed a single representative member of 
each equivalence class. 

We will be examining a number of scenarios of various dimensionalities (that is, differing 
numbers of variables parameterizing p). In all of them, we will seek to find the univariate 
function S scena rio{.v) (our primary computational and theoretical challenge) and the constant 



c 



scenario) 



such that 




(17) 



and 




(18) 
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Given such a pair of volumes, one can immediately calculate the corresponding HS separa- 
bility probability, 

V HS , 

pHS _ sep/ scenario ^ _x 

sep/scenario ~\rHS ' V / 

tot/ scenario 

Let us note that in the full 9-dimensional real and 15-dimensional complex two-qubit 

n 

cases recently studied in [9], it was quite natural to expect that S Tea i(i>) = S rea i(^) (and 
S comp i ex (v) = S comp i ex {^)). But, here, in our /ower-dimensional scenarios, the nullification of 
entries that we employ, breaks symmetry (duality), so we can not realistically expect such a 
reciprocity property to hold, in general. Consequently, we adopt the more general, broader 
formula in (fl3|) as our working formula (TTTj) . 

We now embark upon a series of multifarious lower-dimensional analyses, first for qubit- 
qubit and then qubit-qutrit, qutrit-qutrit and qubit-qubit-qubit systems. These will prove 
useful — as was our original hope — in developing approaches to higher-dimensional analyses, 
presently out of the reach of exact computer analyses. 



II. QUBIT-QUBIT ANALYSES 

To begin, let us make the simple observation that since the partial transposition operation 
on a 4 x 4 density matrix interchanges only the (1,4) and (2,3) entries (and the (4,1) and 
(3,2) entries), any scenario which does not involve at least one of these entries must only 
yield separable states. 



A. Five nullified pairs of off-diagonal entries — 6 scenarios 

1. A- dimensional real case — = f? 

* ^P 1 D 

There are, of course, six ways of nullifying five of the six off-diagonal pairs of entries of p. 
Of these, only two of the six yield any non-separable (entangled) states. In the four trivial 
(fully separable) scenarios, the lower- dimensional counterpart to S rea i(v) was of the form 

^ scenario^) ^scenario 2. 

In one of the two non-trivial scenarios, having the (2,3) and (3,2) pair of entries of p left 
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intact (not nullified), the separability function was 

W") = { (20) 




(It is of interest to note that B u (^, 1) = 2y/p, while in jo], we had conjectured that S rea i{v) 
was proporitional to B u (^, VS)-) 

In the other non-trivial scenario, with the (1,4) and (4,1) pair being the one not nullified, 
the separability function was — in a dual manner (mapping f(y) for v £ [0, 1] into /(£) for 
i/ £ [1, oo]) — equal to 

{2 < v < 1 
(21) 

In both of these scenarios (having c scenar io = 2) for the total (separable and non-separable) 
HS volume, we obtained V* t s = § » 0.0654498 and V™ = ~ 0.0385531. The cor- 
responding HS separability probability for the two non-trivial (dual) scenarios is, then, 
^ » 0.589049. 

16 



5- dimensional complex case — P. 



Now, we allow the single non-nullified pair of symmetrically-located entries to be complex 
in nature (so, obviously we have five variables/parameters — that is, including the three 
diagonal variables — in toto to consider, rather than four). 

Again, we have only the same two scenarios (of the six combinatorially possible) being 
separably non-trivial. Based on the (2,3) and (3,2) pair of entries, the relevant function 
(with the slight change of notation to indicate complex entries) was 




<W^)H (22) 
and, dually, 

{7T < V < 1 
(23) 
I v>l. 

So, the function y/u, which appeared ( (l20j) . (1211 ) in the corresponding scenarios restricted to 
real entries, is replaced by v itself in the complex counterpart. (We note that B u (l, 1) = v.) 
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For both of these complex scenarios, we had V£ t s = ^ and V^f = 3^, for a particularly 
simple HS separability probability of |. 

3. 7- dimensional quaternionic case — P^p = jq 

Here we allow the single pair of non-null off-diagonal entries to be quaternionic in nature 



41 



42| [42, sec. IV]. We found 



v > 1 



2 



and, dually, 



< v < 1 



,T 2 

2^ 2 



1/ > 1. 



(We note that .B„(2, 1) = ^.) For both scenarios, we had V^f = V^ 5 = -^qq , giving 
us pHj~ = jq — which is the smallest separability probability we will report in this entire 
paper. 

So, in our first set of simple (m = 5) scenarios, we observe a decrease in the probabilities 
of separability from the real to the complex to the quaternionic well as a progression 

from yjv to v to v 2 in the functional forms occurring in the corresponding HS separability 
probability functions. 

4- Relevance of Dyson indices 

The exponents of v in the real-complex-quaternionic progression in the immediately pre- 
ceding m = 5 analyses, that is |,1,2 bear an evident elementary relation to the Dyson 



indices 



111 ]. (3 = 1,2,4, corresponding to the Gaussian orthogonal, unitary and symplec- 



tic ensembles 44|. (Further, many of the additional scenarios studied below — also in the 
non-qubit-qubit analyses — will have explicit occurrences in the corresponding separability 
functions of \Jv for real entries and v for complex entries. Of course, use of /x = \J~u as our 
principal variable would give the Dyson series itself, rather than one-half of it.) We note 
that the foundational work of Zyczkowski and Sommers [l|) in computing the HS (separable 
plus nonseparable) volumes itself relies strongly on random matrix theory (in particular, the 
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^aguerre ensemble). Their formula for a certain generalized (Hall) normalization constant 

l[ eq. (4.1)], for instance, contains a dummy variable (3 which equals 1 in the real case and 2 
in the complex case. In their concluding remarks, they write: "these explicit results may be 
applied for estimation of the volume of the set of entangled [emphasis added] states... It is also 

ikely that some of the integrals obtained in this work will be useful in such investigations" 

l|, p. 10125]. 

Of course, random matrix theory is framed in terms of the eigenvalues and eigenvectors 
of random matrices — which do not appear explicitly in the Bloore parameterization — so, it 
is not altogether transparent in what manner one might proceed further to relate the two 
areas. (But for the m = 5 highly sparse density matrices for this set of scenarios, one can 
explicitly transform between the eigenvalues and the Bloore parameters.) 



B. Four nullified pairs of off-diagonal entries — 15 scenarios 



1. 5- dimensional real case — = |; J^z 

Here, there are fifteen possible scenarios, all with V* t s = ^. Six of them are trivial 
(separability probabilities of 1), in which c scenar i is either 7r (scenarios [(1,2), (1,3)], [(1,2), 
(2,4)], [(1,3), (3,4)] and [(2,4), (3,4)]) or 4 (scenarios [(1,2), (3,4)] and [(1,3), (2,4)]). Eight of 
the nine non-trivial scenarios all have — similarly to the 4-dimensional analyses (sec. Ill A II) 
— separability functions S(i>) either of the form, 




Sscenario(L') — \ (^6) 



(for scenarios [(1,2), (2,3)], [(1,3), (2,3)], [(2,3), (2,4)] and [(2,3), (3,4)]) or, dually, 

TV < V < 1 
* V>1 

(for scenarios [(1,2), (1,4)], [(1,3), (1,4)], [(1,4), (2,4)] and [(1,4), (3,4)]). The corresponding 
HS separability probabilities, for all eight of these non-trivial scenarios, are equal to | = 
0.625. This result was, in all the eight cases, computed by taking the the ratio of = ^ 

t yHS _ j± 
lo v tot - 480- 
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In the remaining (ninth) non-trivially entangled case — based on the non-nullified dyad 
[(1,4), (2, 3)] — we have, taking the ratio of V^f = ^ to V^f = a quite different Hilbert- 
Schmidt separability probability of ^ ~ 0.54038. This isolated scenario (with c scenario = 4) 
can also be distinguished from the other eight partially entangled scenarios, in that it is the 
only one for which entanglement occurs for both v < 1 and v > 1. We have 




<%,4),(2,3)](^) = < (28) 



By way of illustration, in this specific case, we have the scenario-specific marginal jacobian 
function, 

rr ( v 07(-3z/ 2 + (^ + 4) + l)log(;/)+3) 

J[(l,4),(2,3)](^ = 30(1/ - l) 5 ' " ' 

2. 6- dimensional mixed (real and complex) case — P^p = i^nfi iH^f! | 

Here, we again nullify all but two of the off-diagonal entries (m = 4) of p, but allow the 
first of the two non-nullified entries to be complex in nature. Making (apparently necessary) 
use of the circular /trigonometric transformation p n = r 2 sin# 2 , p 2 2 = r2 cos9 2 , we were able 
to obtain an interesting variety of exact results. One of these takes the form, 



S IM,M]M = <W,4)]M = ({f ,0 < v < 1} {^v > l) . (30) 

Now, we have V t H J = ^ and V™ = so = ^ « 0.644272. The two dual 

scenarios — having the same three results — are [(1,2), (2,3)] and [(1,3), (2,3)]. 
Additionally, we have an isolated scenario, 



5 [(i;4),(2,3)](^) = {{2^,0 < V < 1} {f,V> 1} 



(31) 



for which, V t * s = ^ and K?/ = Jf^, so Pf* = ±§f « 0.414175. (Note the presence of 
both a/z7 and v in ( 13T1) — apparently related to the mixed [real and complex] nature of this 
scenario (cf. ( 134}) ).) 
Further, 

«%;4),( 2 ,4)]W = 5 [(M) > (3,4) ] (-) = {{f ,0 < v < 1} {g,i/ > 1} , (32) 
the dual scenarios being [(2, 3), (2, 4)] and [(2, 3), (3, 4)]. For all four of these scenarios, 

V HS _ _jf_ ..J T/HS _ _7rf_ pfTS _ 3 _ n oyer 

v tot — 1440 dI1U v sep — 3840' bU sep ~ 8 ~ 
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3. 7- dimensional complex case — P 



sep 



Here, in an m = 4 setting, we nullify four of the six off-diagonal pairs of the 4x4 density 
matrix, allowing the remaining two pairs both to be complex. We have (again observing a 
shift from ^Jv in the real case to v in the complex case) 



S[(1,2),(1~4)]( U ) ~ S[(1~3)A1~A)]( U ' ~ ^[(l^),(274)](^) ~~ "^[(M),(3"4)] 



V) 




Since V t ^ t s = and V^f = y^x), we have P^J? = | = 0.4. We have the same three out- 
comes for the four dual scenarios [(1~2), (2~3)], [(1~3), (2~3)], [(2~3), (2~4)] and [(2~3), (3~4)], 
as well as — rather remarkably — for the (again isolated [cf. (j3"T]) ]) scenario [(1, 4), (2, 3)], hav- 
ing the (somewhat different) separability function (manifesting entanglement for both v < 1 
and v > 1), 



7T 2 U 0<U<1 

^(lAUi^M = { . ( 34 ) 
^ V > 1. 



(However, c scenar io = tt 2 for this isolated scenario, while it equals \ for the other eight.) 
The remaining six (fully separable) scenarios (of the fifteen possible) simply have PfJ^ = 1. 



4- 8 -dimensional mixed (real and quaternionic) case 
We report here that 

c [(i^),(i,4)] = -yg"> c [(i,2),(m)] = 32; ( 35 ) 

where as before the wide tilde notation denotes the quaternionic off-diagonal entry. 

C. Three nullified pairs of off-diagonal entries — 20 scenarios 

1. 6-dimensional real case — P^~ = 2 — |||| 

Here (m — 3), there are twenty possible scenarios — nullifying triads of off-diagonal pairs 
in p. Of these twenty, there are four totally separable scenarios — corresponding to the non- 
nullified triads [(1,2), (1,3), (2,4)], [(1,2), (1,3), (3,4)], [(1,2), (2,4), (3,4)] and [(1,3), (2,4), 
(3,4)] — with c scenario = 2j and V^f = V^f = To proceed further in this 6-dimensional 
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case — in which we began to encounter some computational difficulties — we sought, again, 
to enforce the four nonnegativity conditions ((j2J), ©, (j4]), flS}), but only after setting v to 
specific values, rather than allowing v to vary. We chose the nine values z^ = |,|,|,|, 1, 2, 
3, 4 and 5. Two of the scenarios (with the triads [(1,2), (2,3), (3,4)] and [(1,3), (2,3), (2,4)]) 
could, then, be seen to fit unequivocally into our earlier observed predominant pattern, 
having the piecewise separability function, 

^ < < 1 

«%,2),(2,3),(3,4)]0) = «%,3),(2,3),(2,4)](^) = { (36) 

\ U>1. 

We, then, computed for these two scenarios that V t % s = gfgo ~ 0.00538303 and (again mak- 
ing use of the transformation pn = r 2 sin 9 2 , p 22 = r 2 cos 9 2 ) that V^f = 2 ^gfgQ — j^ll 2 
0.00358207. This gives us P s f/ = 2 - ||f « 0.665437. For two dual dyads, we have the 
same volumes and separability probability and, now, the piecewise separability function, 

T - v - 1 

<%,2),(1,4),(3,4)] (f) = «5[(i )3 ) ) (i ) 4) ) (2,4)](i / ) = "J (37) 

We have not, to this point, been able to explicitly and succinctly characterize the functions 
<S scenario^) for non-trivial fully real m = 3 scenarios other than the dual pair ( ([561) . (I3T|) ). 

In all the separably non-trivial scenarios so far presented and discussed, we have had 
the relationship iS scenario (l) = c scenario . However, in our present m = 3 setting (three pairs 
of nullified off-diagonal entries), we have situations in which S scenar i (l) < c scenar i . The 

2 

values of c scenar i in the sixteen non-trivial fully real m = 3 scenarios are either ^- ?a 4.9348 
(twelve occurrences) or ^ w 4.18879 (four occurrences— [(1,2), (1,3), (1,4)], [(1,2), (2,3), 
(2,4)], [(1,3), (2,3), (3,4)] and [(1,4), (2,4), (3,4)]). In all four of the latter (^) occurrences, 
though, we have the inequality, 

1 47T 
S scen ario{l) = ^ ( 12 + + ^) « 3 ' 8281 < "J ~ 4.18879, (38) 

as well as a parallel inequality for four of the twelve former (^-) cases. The implication 
of these inequalities for those eight scenarios is that at v = 1 (the value associated with 
the fully mixed [separable] classical state), that is, when puPu = P22P33, there do exist 
non-separable states. 
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2. 7- dimensional mixed (one complex and two real) case — P^p = 

Here, in an m = 3 setting, we take the first entry of the non-nullified triad 
to be complex and the other two real. Of the twenty possible scenarios, four 
- [(1~2), (1, 3), (1, 4)], [(1,2), (2, 3), (2, 4)], [(1~3), (2, 3), (3, 4)] and [(1~4), (2,4), (3,4)]— had 

2 

^scenario = \ ~ 4.9348 and these four all had the same (lesser) value of 

56 7T 2 

SscenarM = — + — « 4.54148. (39) 

There were seven scenarios with c scenar io = ^ ~ 5.58505. Three 

of them— [(1^2), (1,3), (2,4)], [(1,2), (1,4), (2,3)] and [(1,3), (1,4), (2,3)] 

—had S scenar i (l) = (manifesting equality), while four— 

[(1,2), (1, 3), (2, 3)], [(^2), (1,4), (2, 4)], [(1,3), (1,4), (3, 4)] and [(2~3), (2, 4), (3, 4)]— had 
the result (!39|) (manifesting inequality). 

The remaining nine of the twenty scenarios all had c scenar i = S scenaT i {l) = ^j- ~ 6.57974. 
For one of them, we obtained 

%;2),( 2 ,3),(3,4)]M= 2i> <^<i (40) 

with associated values of V t ^ t s = t^q, V^ p s = anc ^ ^sep = I§ ~ 0.6875. A dual 

scenario to this one that we were able to find was [(1, 2), (1, 4), (3, 4)]. The separability 
functions — and, hence, separability probabilities — for the other eighteen scenarios, however, 
are unknown to us at present. 



3. 8-dimensional mixed (two complex and one real) case 
Our sole result in this category is 

c [(i;2),(i; 3 ),(i,4)] - 15 • l 4i J 

4- 9-dimensional complex case 

Now, we have three off-diagonal complex entries, requiring six parameters for their spec- 
ification. This is about the limit in the number of free off-diagonal parameters for which we 
might hopefully be able to determine associated separability functions. 
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As initial findings, we obtained 

5 [(i;4) i ( 2 ;3),(2;4)]( 1 ) = c [(i;4),(2; 3 ),(2;4)] = y> ( 42 ) 

and also for scenarios [(1~2), (1~4), (3~4)], [(1~3), (1~4), (2~4)] and [(1~4), (2~3), (3~4)], while 

7T 3 

c [(i^),(i;3),(M)] = c [(i;3),(2; 3 ),(3;4)] = c [(i;4),( 2 ;4),( 3 ;4)] = y • ( 43 ) 

D. Two or fewer nullified pairs of off-diagonal entries 

1. 7- dimensional real case 

The [(1,2), (1,3), (2,4), (3,4)] scenario is the only fully separable one of the fifteen possible 
(m = 2). For all the other fourteen non-trivial scenarios, there are non-separable states both 
for v < 1 and v > 1. For all fifteen scenarios, we have c scenar io = ~ 6.57974. Otherwise, 
we have not so far been able to extend the analyses above to this m = 2 fully real case (and 
a fortiori the m = 1 fully real case), even to determine specific values of S scenario {l). 

2. 8 -dimensional real case 

Here we have c scenar i = ~ 8.77298 for all the six possible (separably non-trivial) 
scenarios (m = 1). Let us note that this is, in terms of preceding values of these constants 
(for the successively lower-dimensional fully real scenarios), ^ = |(^f-), while ^ = f (^r)- 
Also, = |(^p), the further relevance of which will be apparent in relation to our 
discussion of the full 9-dimensional real scenario (sec. IIXI ). 



III. QUBIT-QUTRIT ANALYSES 

The cancellation property, we exploited above, of the Bloore parameterization — by which 
the determinant and principal minors of density matrices can be factored into products 
of (nonnegative) diagonal entries and terms just involving off-diagonal parameters (%)— 
clearly extends to n x n density matrices. It initially appeared to us that the advantage of 
the parameterization in studying the two-qubit HS separability probability question would 
diminish if one were to examine the two-qubit separability problem for other (possibly mono- 
tone) metrics than the HS one (cf. [8|), or even the qnbit- qutrit HS separability probability 
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question. But upon some further analysis, we have found that the nonnegativity condition 
for the determinant of the partial transpose of a real 6x6 (qubit-qutrit) density matrix (cf. 
d2J)) can be expressed in terms of the corresponding zy's and two ratio variables (thus, not 
requiring the five independent diagonal variables individually), 

v\ = , u 2 = , (44) 

P22P44 P33P55 

rather than simply one (y) as in the 4x4 case. (We compute the qubit-qutrit partial 
transpose by transposing in place the four 3x3 blocks of p, rather than — as we might 
alternatively have done — the nine 2x2 blocks.) 



A. Fourteen nullified pairs of off-diagonal entries — 15 scenarios 

1. 6- dimensional real case — Pj%£ = ff 

To begin our examination of the qubit-qutrit case, we study the (m = 14) scenarios, in 
which only a single pair of real entries is left intact and all other off-diagonal pairs of the 
6x6 density matrix are nullified. (We not only require that the determinant of the partial 
transpose of p be nonnegative for separability to hold — as suffices in the qubit-qubit case, 



given that p itself is a density matrix 



45 



46] — but also, per the Sylvester criterion, a nested 



series of principal leading minors of p.) We have six separably non-trivial scenarios. (For all 
of them, V£ 3 = ^0 

Firstly, we have the separability function, 




The dual scenario to this is [(2,4)]. Further, 



2 



'[(1.6)] 

with the dual scenario here being [(3,4)]. Finally 



VlV 2 < 1 
V\V 2 > 1, 



2 v 2 < 1 



-7= v 2 > 1, 
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having the dual [(3,5)]. 

The remaining nine possible scenarios — the same as their complex counterparts in the 
immediate next analysis — are all fully separable in character. 

We have found that V^f = for the six non-trivially separable scenarios here, so 
P^s — || pa 0.589049, as in the qubit-qubit analogous case (sec. Ill A 1H . 

2. 7- dimensional complex case — -P^f = \ 

Now, we allow the single non-nullified pair of off-diagonal entries to be complex in nature 
(the two paired entries, of course, being complex conjugates of one another). {V^f = 
for this series of fifteen scenarios.) Then, we have (its dual being [(2,4)]) 



S^M = < - (48) 




Further, we have (with the dual [(3,4)]) 



, 7T V\Vn<\ 

S^ )] (u 1 ,u 2 )={ ' (49) 



- 21 - V\Vi > 1 



and (its dual being [(3,5)]), 



, 7T Po < 1 

s ^ = - (so) 

For all six of these scenarios, V sep = j^d, so P^g = |. 

B. Thirteen nullified pairs of off-diagonal entries — 105 scenarios 

i. 7- dimensional real case — = |; || ; j4r 

Continuing along similar lines (m = 13), we have 105 combinatorially distinct possible 
scenarios. Among the separably non-trivial scenarios, we have 



6x6 6x6 J 71 Vl ~ 1 

5 Kl,2),(l,5)]( I/ l) = <5 [ 6 (l,4),(l,5)]( Z/ l) = "! ( 51 ) 
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(duals being [(1,2),(2,4)] and [(1,4),(2,4)]). We computed V* s = ^,V S % S = ^ so 
— I — 0.625 for these scenarios. 



Also, 



7T V\l>2 < 1 

5 [(1,3),(1,6)]( Z/ 1' ^2) = « 5 p,4),(l,6)]( Z/ l' ^2) = <j (52) 

^77^ ^2 > 1. 



We, then, have V t » s = V s % s = so Pf* = A = 0.3125. 

Additionally, 



,S U1A).(2.6)]( U ^ ~ { ( 53 ) 



4 i/ 2 < 1 



1(1,4),(2,6)] 



4= ^ > 1. 



For this scenario, we have V^ 5 = V£/ = so P £jf = 1 ~ 0.294524. 

Further still, 

«),( 2 ,4),K)= P ^ 
[4^ Z/!<1. 

For this scenario, we have V£ t s = V s % s = ^5, so Pf* = 3^ « 0.54038. 
Further, 

{7T f 2 < 1 
2 ( cos (\A-^) + ^ 1 ) i/2>1 - 

The separability function for [(1,3), (1,5)] is obtained from this one by replacing u 2 by v 1 . 
Also, 

5 ?a, 6 2),(3^)](^^)= ( 56 ) 

4./1--* * V Z. z/ lZ / 2 >l. 



The separability function for [(1,2), (3, 5)] can be obtained from this one by setting v\ = 1. 



1057T 



2. 8- dimensional mixed (real and complex) case — P^f = 51 , 
Further, we have (with V t ^ t s = for all scenarios), 



■ 4? 1/1 > 1 

76x6 /,, \ _ 06x6 /,, \ _ ) 3 



^w*>=«w*>=w „ , (57) 

X < ^ < 1. 
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sep . f unct . 




FIG. 1: Plot of the separability function S 6 *-? , (v x ) 

[(2,3), (2,4)] V L/ 



Since V£? 



sep 393216 

Further, 



we have P£f 



_ 1057T 

sep 512 



<-6x6 / > 

6 [(0),(2,4)]^ 



0.644272 for both these scenarios. 



^ < v x < 1 



and 



[(l,3),(3,4)]^ l! 2) [(1, 

Additionally, 



6x6 / 
[(M),(3,4)] ( 
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3 



^1^2 > 1 



^^1^2 < ^1/2 < 1. 



(58) 



(59) 



and 



5 [ 6 (S),(3,4)](^ 2 ) 



in 

3 



^1^2 > 1 



\^^JV\V2 (3 — ^1^2) < Z/i^ 2 < 1 



iZL 

3 



|7T V /77I^ 



Z/iZ/2 = 1 

< U X U 2 < 1 



+ u x u 2 >l. 



We have also obtained the separability function (Fig. [T]) 



(60) 



(61) 



S bxb lux) 

[(2,3),(2,4)] V X > 



3 



V\ > 1 



\tt (3 - vi) y/u{ < u x < 1. 



(62) 



Of the 105 possible scenarios, sixty had S! 



6x6 Q ^ 

scenariox 1 I 



45 

scenario 3 



thirty-three had 



3 



6x6 ■ (1 1) = C 

scenano\ L i L ) ^scenario 

2tt. 



2tc, and twelve (for example, [(3,4), (5,6)]) had 5" 



6x6 /i 
scenario V ) / 
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3. 9- dimensional complex case — P^ s — -■ - 



sep 3) 5 



We have obtained the results 



_2, 



- v x > 1. 

Since V t ^ t s = 36 ^ 80 and Vj^f = 90 jl 00 , we have here = | = 0.4. We have the same 
three outcomes also based on the separability function, 

I TT 1/1 1/9 > 1 

S 6x6 - (^,i/ 2 ) = ^ 2 (64) 

[(1,4),(3,4)] V ^ ^ \ 



1^2 



7T i^i ^2 V\l>2 < 1. 



Further, 



and 



, 2 ^^-7T 2 jy jy > 1 

*sw^>- ,r „ <, <65> 

±7rVii/ 2 < viv 2 < 1. 



< I/x < 1 



* - 5T "1 > 

For both of these last two scenarios, we have V^ 5 = 362 g 80 and V^£f = 362 g 80 , leading to 
Flip = I ~ 0.33333. Also, we have these same three outcomes based on the separability 
function, 

, , UiUo > 1 

^w,~ 4 >i^)= (67) 



[(2,3),(3,4)]^ 



\-K 2 ViV 2 (2 - ^1^2) < z/iZ/ 2 < 1- 



2 

Of the 105 possible scenarios — in complete parallel to those in the immediately preceding 
section— sixty had S^ e l ario (l, 1) = c scenario = ^ , thirty-three had Sl^ nario {l, 1) = c scenario = 
7T 2 , and twelve (for example, [(3,4), (5,6)]) had S^ e l ario (l, 1) = ^ < Cscenario = k 2 . 

Our results in this (9-dimensional) section and the (8-dimensional) one immediately pre- 
ceding it are still incomplete with respect to various scenario-specific separability functions 
and, thus, the associated HS separability properties. 
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C. Twelve nullified pairs of off-diagonal entries — 455 scenarios 



1. 8- dimensional real case 



Now, we allow three of the off-diagonal pairs of entries to be non-zero, but also require 
them to be simply real. We found the separability function 



<5[(l,2),(l,3),(3,4)]( I/ l, I/ 2) — <5[(l,2),(l,4),(3,4)]( Z/ 'l, l/ 2) 



(68) 



■J.7T 

3 



v 2 A v\ > 



§7T ^3 SBC" 1 (y/UiU2) + 4:y/Vih> 2 + 

Also, we have 



\/v\V2 — 1 



v x > A j- > v 2 A v 2 > 
4\M^2 - 1 ) ^i>0Ai<f 2 . 



<%,3),(l,4),(2,4)](^l) = < 



4-7T 

3 

47r v ^75" 



7T (4 ! / 1 3/2 + (3 sin- 1 (^/l- -4 v / ^T r T) ^i+v 7 ^! 3 ! 



3fi 



^ = 1 

< v x < 1 
1/1 > I. 



(69) 



We note, importantly, that in all the qubit-qutrit scenarios in which v\ and v 2 have both 
appeared in the (naively, bivariate) separability function, it has been in the 'product form 
v\v 2 (cf. sec. IX Ap . 



IV. QUTRIT-QUTRIT ANALYSES 

In the qubit-qubit (4x4 density matrix) case, we were able to express the condition (JSJ) 
that the determinant of the partial transpose of p be nonnegative in terms of one supplemen- 
tary variable (is), given by (jSJ), rather than three independent diagonal entries. Similarly, in 
the qubit-qutrit (6x6 density matrix) case, we could employ two supplementary variables 
(^lj v 2)i given by (1441) . rather than five independent diagonal entries. 

For the qutrit-qutrit (9x9 density matrix) case, rather than eight independent diagonal 
entries, we found that one can employ the four supplementary variables, 

PllP55 P22P66 P44P88 P55P99 , nn \ 

v x = ; v 2 = ; u 3 = ; i/ 4 = . (70) 

P22P44 P33P55 P55P77 P66P88 
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A. Thirty-five nullified pairs of off-diagonal entries — 36 scenarios 



1. 10- dimensional complex case — Ppp? = |i | 

Here, we nullify all but one of the thirty-six pairs of off-diagonal entries of the 9x9 
density matrix p. We allow this solitary pair to be composed of complex conjugates. Since 
the Peres-Horodecki positive partial transposition (PPT) criterion is not sufficient to ensure 
separability we accordingly modify our notation. 

Our first result is 



, 7T Vi < 1 



(a dual scenario being [(2,4)]). We have V^ 5 = gg^soo, V^ 5 T = 108 ^ 400 , so P^/ T = ±. 
The same three outcomes are obtained based on the PPT function 



, 7T Z/iZ/9 < 1 

. " (72) 

^1^2 > 1- 



On the other hand, we have V^ 5 = 362 g 800 , ^p^r = 21772800 ' anc ^ -^ppt = | based on the 
PPT function 

= """" < 73 > 

Of the thirty-six combinatorially possible scenarios, thirteen had Ppp T = |, while four 
had PppT = k, and the remaining nineteen were fully separable in nature. 




B. Thirty- four nullified pairs of off-diagonal entries — 630 scenarios 

1. 12- dimensional complex case — PppT = §! 35 

Since the number of combinatorially possible scenarios was so large, we randomly gener- 
ated scenarios to examine. 
Firstly, we found 

f7T 2 1/2 > 1 

. „" , (74) 
ir 2 is 2 < vi < 1. 
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For this scenario, we had V t * s = 479 ^ W00 , V^ T = 1437 ^ 4 8oo ^ S ivin g us P ppt = §■ 
Also, we found 

v 2 v 4 < 1 

„ „ 1 x <75 » 

For this scenario, we had V t * s = 479 ^ imo , V^ T = 2052 ll 4000 , giving us P^ T = ^ « 0.23333. 



V. QUBIT-QUBIT-QUBIT ANALYSES, I 

For initial relative simplicity, let us regard an 8 x 8 density matrix p as a bipartite system, 
a composite of a four-level system and a two-level system. Then, we can compute the partial 
transposition of p, transposing in place its four 4x4 blocks. The nonnegativity of this partial 
transpose can be expressed using just three ratio variables, 

PllP66 P22P77 P33P88 ,„ c s 

v\ = ; vi = ; vz = , (76) 

P22P55 P33P66 P44P77 

rather than seven independent diagonal entries. 



A. Twenty-seven nullified pairs of off-diagonal entries — 28 scenarios 

1. 9- dimensional complex case — Ppp^ 
We have the PPT function 




v x < 1 

(77) 



(Scenario [(2,5)] was dual to this one.) For this scenario, V t ^ t s = 362 ^ 80 , Vpp T = 108 g 640 ; 
yielding Ppp T = |. There were twelve scenarios, in toto, with precisely these three outcomes. 
The other sixteen were all fully separable in nature. 



B. Twenty-six nullified pairs of off-diagonal entries — 378 scenarios 

1. 11- dimensional complex case — Pppx = |i § 



Again, because of the large number of possible scenarios, we chose them randomly for 
inspection. 

28 



Firstly, we obtained 



\^M [Ul ^ 2) -\ 2 nA 1 a n 

tx 1 v\v 2 V\ > A > v 2 A v 2 > 0. 



(78) 



(Of course, the symbols "A" and "V", used by Mathematica in its output, denote the logical 
connectives "and" (conjunction) and "or" (intersection) of propositions.) For this scenario, 

we had V t ^ t s = 399 i 6800 , Vppr = 119750400 > g ivin g us Pppr = §■ 
Also, 



7T 2 


V\ > 1 A V3 > 1 




< f 1 < 1 A 1/3 > 1 


7T 2 Z/ 3 


y 1 >lA0<j/ 3 <l 




< v x < 1 A < v 3 < 1 



For this scenario, we had V^ t s = 3g9 f 6S00 , V^ T = 35g ^ 1200 , giving us P^ T = §. 
We also found the PPT function 



2 



V 3 = 1 

< 2/3 < 1 



7T COS 



-1 



1 - i 

^3 



sm 



^3 + 7T 2 - ^3>1 



(79) 



(80) 



VI. QUBIT-QUBIT-QUBIT ANALYSES. II 



Here we regard the 8x8 density matrix as a tripartite composite of three two-level 
systems, and compute the partial transpose by transposing in place the eight 2x2 blocks of 
p. (For symmetric states of three qubits, positivity of the partial transpose is sufficient to 
ensure separability [47|, |48j.) Again the nonnegativity of the determinant could be expressed 
using three (different) ratio variables, 



pupu , 
P22P33 



PuP55 , 
P33P66 



^3 



P55P88 

PmPn 



(81) 
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A. Twenty-seven nullified pairs of off-diagonal entries — 28 scenarios 



9- dimensional complex case — Pppp = 3 



There were, again, twelve of twenty-eight scenarios with non-trivial separability proper- 



ties, all with V t 



HS 

tot 



362880 ' 



V PPT = Tns!fi4n> y ieldin g P PPT = \- 0ne of these was 



7T V\ < 1 



^ V X > 1. 



(82) 



11- dimensional complex case — Pppt = M> 3 



We obtained the PPT function 



7Y 



^ = u 3 A v x > A v 2 > 

v 2 > A ((i/! = A i/ 3 > 0) V (u 3 = A v x > 0)) 
1/1 > A v 2 > A ^ < 1/3 



vi 



n 2 _ ^3 i/ x > A 1/2 > A ^ > 1/3 A 1/3 > 0. 



VI 



For this we had V t g s = og^o, V™ t = 



17tt 2 



239500800 



Additionally, 



[(M),(W U 



7T 



1/1 < 1 

V 1 >1. 



Here, we had V^ 5 



y #PT = 1197^00 ' giving us P^/ T 



tot 39916800 ' 

Another PPT function we were able to find was 



(83) 



, giving us P^ T = g « 0.283333. 



(84) 



5 [(S)W^ 3) 



7T 

2 



v 2 > A (1/3 = i>2 V (z/ 2 > ^3 A z/ 2 < 2^3) V {y 2 > 2z/ 3 A v 3 > 



(85) 
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VII. APPROXIMATE APPROACHES TO 9-DIMENSIONAL REAL QUBIT- 
QUBIT SCENARIO 



As we have earlier noted, it appears that the simultaneous computational enforcement of 
the four conditions ((j2J), Q, (jlj), (jSJ)) that would yield us the 9-dimensional volume of the 
separable real two-qubit states appears presently highly intractable. But if we replace (jSj) 
by less strong conditions on the nonnegativity of the partial transpose (p T ), we can achieve 
some form of approximation to the desired results. So, replacing (jSJ) by the requirement 
(derived from a 2 x 2 principal minor of p T ) that 

1 - uz 2 u > 0, (86) 

we obtain the approximate separability function 

IW 0<i/<1 

(In the anises in this section, we utilize the integration litnits on the V . fl e qS . (3)-(5)] 
yielded by the cylinrical decomposition algorithm [CAD], to reduce the dimensionalities of 
our constrained integrations.) This yields an upper bound on the separability probability 
of the real 9-dimensional qubit-qubit states of \ + 1 ^ 2 ~ 0.88427. We obtain the same 
probability if we employ instead of (186)) the requirement 

v - z 2 23 > 0, (88) 

which yields the dual function to (187)) . namely, 



512tH 



V > 1 



S%T(") = i 27 (89) 
^ 2 (3-z/)v^ 0<i/<1. 

(The left-hand sides of ( 1861) and (188)) are the only two of the six 2 x 2 principal minors 
of p T that are non-trivially distinct — apart from cancellable nonnegative factors — from the 
corresponding minors of p itself.) The non-constant functional form in the second line of 
( 1891) will emerge again, importantly, in ( 193)) . 

If we form a "quasi-separability" function over v G [0, oo] by piecing together the 
non-constant segments of (187)) and (1891 . we can infer — using a simple symmetry, dual- 
ity argument — an improved (lowered) upper bound on the HS separability probability of 
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- 0.76854. We can also reach such a result by noting that the two constraints (!86|) 
and (j88p are independent (involve different variables), so we should just be able to multiply 
the corresponding functions (and then scale them by the corresponding c scenar i = 51 |J 2 ) 



VIII. ALTERNATIVE USE OF BLOCH PARAMETERIZATION 

We may say, in partial summary that we have been able to obtain certain exact two-qubit 
HS separability probabilities in dimensions seven or less, making use of the advantageous 
Bloore parameterization [h]] , but not yet in dimensions greater than seven. This, however, is 
considerably greater than simply the three dimensions (parameters) we were able to achieve 



17J in a somewhat comparable study based on the generalized Block representation param 



eterization 



50] . In [17|] — extending an approach of Jakobczyk and Siennicki 50| — we 
primarily studied two-dimensional sections of a set of generalized Bloch vectors correspond- 
ing tonxn density matrices, for n = 4, 6, 8, 9 and 10. For n > 4, by far the most frequently 
recorded HS separability [or positive partial transpose (PPT) for n > 6] probability was 
f ~ 0.785398. A very wide range of exact HS separability and PPT probabilities were 
tabulated. 

Immediately below is just one of many matrix tables (this one being numbered (5)) 
presented in [17J (which due to its copious results has been left simply as a preprint, rather 
than submitted directly to a journal). This table gives the HS separability probabilities 
for the qubit-qutrit case. In the first column are given the identifying numbers of a pair 
of generalized Gell-mann matrices (generators of SU(6)). In the second column of fl90l) are 
shown the number of distinct unordered pairs of SU (6) generators which share the same 
total (separable and nonseparable) HS volume, as well as the same separable HS volume, 
and consequently, identical HS separability probabilities. The third column gives us these 
HS total volumes, the fourth column, the HS separability probabilities and the last (fifth) 
column, numerical approximations to the exact probabilities (which, of course, we see — 
being probabilities — do not exceed the value 1). (The HS separable volumes too can be 
deduced from the total volume and the separability probability.) 
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{1,13} 
{3,11} 4 
{3,13} 4 
{3,25} 4 
{8,13} 4 
{8,25} 4 
{11,15} 4 
{11,24} 2 
{13,24} 2 
{13,35} 4 
{15,16} 4 
{16,24} 2 
{20,24} 2 
{24,25} 2 
{24,27} 2 ^ 
{25,35} 4 
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V3tt 
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1 5^5 
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80^ 
V / 5+3csc- 1 ( 




3tt 



0.785398 
0.707107 
0.833333 
0.883883 
0.577350 
0.816497 
0.746830 
0.863648 
0.979236 
0.886838 
0.879838 
0.549815 
0.685627 
0.821115 
0.903076 
0.504975 
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It might be of interest to address separability problems that appear to be computationally 
intractable in the generalized Bloch representation by transforming them into the Bloore 
parameterization. 



IX. FULL REAL AND COMPLEX TWO-QUBIT SEPARABILITY PROBABILITY 
CONJECTURES 

The qubit-qubit results above (sec. [Ill) motivated us to reexamine previously obtained 



results (cf. 
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eqs. (12), (13)]) and we would like to make the following observations 
pertaining to the full 9-dimensional real and 15-dimensional HS separability probability 
issue. We have the exact results in these two cases that 

J Jreal(v) = 2 J J real {p) = « §.60561 • 10~ 6 (91) 

and 

JcompUv) = 2 J J comp leM = 10Q9008Q00 « 9 - 91072 • 10_1 °- ( 92 ) 
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Now, to obtain the corresponding total (separable plus nonseparable) HS volumes computed 
by Zyezkowski and Sommers that is, ~ 0.0016106 and 8513 ^ 6 0500 ~ 1.12925 ■ 10~ 6 , 

one must multiply (ED) and ([92]) by the factors of C reai = = ^ « 187.157 and 

C CO mpiex = = ^f - ~ 1139.42, respectively. 

To most effectively compare these previously-reported results with those derived above 
in this paper, one needs to multiply C rea i and S rea i(is), by 2~ 4 = ^ and in the complex case 
by 2~ 7 = y^g. Doing so, for example, would adjust C rea i to equal c rea i = ^§f^ ~ 11.6973, 
which we note, in line with our previous series of calculations [sec. ITTTJ2] is equal to f(^). 
(Andai [s| also computed the same volumes — up to a normalization factor — as Zyezkowski 
and Sommers [lj].) Now, our estimates from j^] are that S rea i(l) = 114.62351 < C rea i and 
ScompiexiX) = 387.50809 < C comp i ex . These results would appear — as remarked above — to 
be a reflection of the phenomena that there are non-separable states for both the 9- and 
15-dimensional scenarios at v — 1 (the locus of the fully mixed, classical state). 

Alternatively, the results in sec. Ill A 31 and further throughout the paper, in which we 
find a relation between separability functions and the Dyson indices ((3 = 1, 2, 4) of random 
matrix theory — including the frequent occurrence of y/v in a real scenario and v in the 
corresponding complex scenario — strongly suggest that in the full (m = 0) 9-dimensional 
real and 15-dimensional complex cases scenarios, the separability function for the complex 
case might simply be proportional to the square of the separability function for the real case 
(and, in the quaternionic case j^], to the fourth power of that function). 

Following such a line of thought, we were led to reexamine the numerical analyses re- 
ported in |9J] , in which we had formulated our beta function ansatze. In Fig. [2] we show 
the previously-obtained numerical estimates of S^iiy) and S comp i ex (iy), now both scaled 
("regularized") to equal 1 at v — 1 ., along with the similarly regularized form (termed the 
"incomplete beta function ratio" [40] or, alternatively, the "regularized incomplete beta 
function" ) 

J w (i/, i, 2) = i(3- V)yfr (93) 

of the incomplete beta function, B v {y^i) = |(3 — v)\fv) and l v {v, |,2) 2 . (Let us make 
the important observation here that the functional form ( |93|) has — up to proportionality — 
already occurred [although we did not immediately perceive then its beta function expres- 
sion] in certain previous exact qubit-qubit analyses (189]) (see also (1601) (|62l . but also (1671) . 
for its occurrence in the qubit-qutrit context). 
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FIG. 2: The most subordinate of the three curves is actually two virtually indistinguishable curves: 
(1) our normalized (previously obtained [9|]) numerical estimate of S cornp i ex (u); and (2) the (ex- 
traordinarily well- fitting) square of the incomplete beta function ratio I v (y, ^,2) = i(3 — v)^Jv. 
The intermediate curve is the normalized previously obtained numerical estimate of S rea i{v) and 
the rather well-fitting (most dominant) curve is Iv(v, ^,2). 

Fig. [2] does reveal an extraordinarily good fit between the normalized numerical estimates 
of Scompiex(v) and l v {y, §, 2) 2 , while I u (u,^,2) itself provides a close fit to the normalized 
numerical estimates of S rea i(u). (Note that I u (v,^,2) does contain a factor of ^Jv and 
I v (y, ~, 2) 2 , obviously a factor of u, much in line with the more elementary lower- dimensional 
real-complex examples studied in sec. [Til Further, as an exercise, we sought that value of x for 
which the function I v (i>,^,x) would when employed in our basic paradigm here, as in Fig. [2, 
jointly minimize the sum of a certain least-squares fit to the normalized numerical estimates 
of S rea i(v) and S comp i ex (u). Our Mathematica program produced the answer x = 1.88487, 
being somewhat intermediate in value between ~ 1.732 and 2, the exact candidate values 
we have considered, which both fit the numerical results of rather well.) So, it would seem 
appropriate to revise the two central ansatze put forth in |9j to account for these interesting 



newly-observed phenomena inherent in the results already reported in 
We can now exactly perform the requisite integrations (cf. (flTl) . (JT8 



9] 



r ill 

2/ J rea i(u)IJu, ~,2)dv = = — , (94) 

J v ; v '2' ' 151200 2 5 -3 3 -5 2 -7' V ; 

1 1 71 71 

2 / JcamplexWM ^ 2 ) 2 ^ = gggg^QOO = V ■ 3 4 • 5 3 • 7 2 • ll 2 ■ 13 ' (95) 







f 1 1 5989 

2 ^ J qU atemionicWu{^ ^ 2 ) 4(iu = 358347086242825680000 (%) 
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53- 113 

2 7 • 3 4 - 5 4 - 7 2 • ll 2 • 13 2 • 17 2 • 19 2 • 23 2 ' 
The three marginal univariate jacobian functions above are obtained by transforming the 

3/3 

jacobian for the Bloore parameterization — (nf =1 pa)~2 , (3 = 1,2,4 — to the v variable and 
integrating over the two remaining independent diagonal entries of p. 

So, assuming the validity of our modified beta function ansatze for the real and com- 
plex separability functions, all we still lack for obtaining the Hilbert-Schmidt separable 
volumes/probabilities of the 9-dimensional real and 15-dimensional complex qubit-qubit 
systems themselves are the appropriate (presumptively, exact in nature) scaling constants 
(on the order of 114.61 and 387.467 jsj]) by which to multiply the results of ( 194)) and ( 1951) . 
(We will presume — in light of the numerous analyses reported earlier — that such scaling 
constants are exact in nature, being of the form where k are natural numbers. We 
search over the spaces of possibilities to find choices that accord with our previously-obtained 
numerical results for the HS separability probabilities.) 



A. Real Two-Qubit Case 

If we employ ~ 114.599 as the scaling constant in the real case — giving us a very 
good fit to the numerical estimate of S rea i(l) ~ 114.61 — we obtain an HS separable volume 

of 



-, 2 g 52Q and an HS separability probability of A ~ 0.470588. (Using the numerical results 



9j, we were able to obtain an estimate of this probability as close as 0.46968 by replacing 
the jacobian function (TT5T) by a sixth-order Taylor series approximation of it around v — ^ . 
Providing inferior fits to 114.61, but still of possible interest, would be choices of scaling 
constants ~ 113.644 and 2pl « 115.448. These would lead to HS real separability 
probabilities of ^ ~ 0.466667 and ^ ~ 0.474074 — with the first of these two seeming much 



15 v.^ww 135 
more consistent with the numerics of 



9j than the second. 



By the ("twofold") theorem of Szarek, Bengtsson and Zyczkowski |20j (cf. [5l|] ) — 
formalizing results in [18] — the HS separable volume of the generically rank-3 real qubit- 
qubit states would — adopting as the appropriate scaling constant — be and the 
HS separability probability, ^ w 0.235294. (The HS area-volume ratio for the 9-dimensional 
real two-qubit states is 18^ ~ 31.1769 [l|, eq. (7.9)], while the analogous ratio restricted 
to the separable subset is one-half as large, that is, 9v / 3 ~ 15.5885, indicating the more 
hyperspherical-like shape of the separable subset). 
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B. Complex Two-Qubit Case 



One simple candidate for the scaling constant in the complex case is 2|i « 384.566. 
This would yield an HS separability probability of ||| ~ 0.242045. But considerably more 
attractive (certainly, in part, due to its interesting consonance with the real results just 
advanced, and also the presence of 256 = 2 s and 639 = 9-71, with the 71 in (1951) . thus, 
being cancelled), it seems is ^§§- ~ 385.157 (slightly closer also to our estimate of 387.467 
from |9j). This choice would yield an HS separable volume of ^mlims ~ 2.73758 • 10~ 7 , 
and separability probability of Jj ~ 0.242424, very close to our previous numerically-derived 



estimate of 0.242575 (implicitly given in 18|, between eqs. (41) and (42)]) (and only slightly 
more than one-half of yy). (Let us also indicate that the HS area- volume ratio for the 15- 
dimensional complex two-qubit states is 30-^3 ~ 51.9615 fl, eq. (6.5)], while the analogous 
ratio restricted to the separable subset is again one- half, that is 15^3 « 25.9808, the lesser 
value indicating the more hyperspherical-like shape of the separable subset). 

In the real and complex analyses just conducted, we have tacitly assumed — as we will 
also do in the succeeding, remaining ones — that the appropriate scaling constants should 
be of the form —j-, where, in addition to i and j being natural numbers, k is identical to 
the power of n occurring in the Zyczkowski-Sommers/Andai formulas for the corresponding 
total volumes. Doing so, at least seems plausible, in light of our numerous lower- dimensional 
analyses above. 

The simplicity of two-qubit complex and real HS separability probabilities, J| and ^, 
apparently stemming from the use of the modified beta function ansatze, now leads us to 
examine if we can generate somewhat parallel HS separability probabilities for the 15- and 
35-dimensional real and complex qubit- qutrit cases. 



X. FULL REAL AND COMPLEX QUBIT-QUTRIT SEPARABILITY PROBABIL- 
ITY CONJECTURES 

A. Real Qubit-Qutrit Case 

In Fig. [3] we show an interpolated estimate (with parameter values restricted to the unit 
square) of the real qubit-qutrit separability function. (Auxiliary analyses give very strong 
evidence, as certainly seems plausible, that this function is symmetric under the interchange 
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FIG. 3: Interpolated estimate over the unit square of the real qubit-qutrit separability function 
Srealiyii v t), based on 785,000 randomly generated 6x6 real density matrices 

of v\ and u 2 .) An immediate conjecture, suggested by our various earlier qubit-qutrit results 
(sec. II III and f l60|) ) is that this (naively, bivariate) function is actually univariate in nature, 
and satisfies the proportionality relation (cf. (!93l) ) 

Sreal(vi, = Sreal(v) °C I^V, ~> 2) = ^(3 - rj)<ffj. (97) 

Here 77 = v\v 2 = pi ^f 6 ^ ; being independent of p 22 and P55 (given the definitions of v\ and v 2 
in eq. (jSJ) above). (If we had chosen to compute the partial transpose of p by transposing 
in place its nine 2x2 blocks, rather than its four 3x3 blocks, then presumably the same 
essential phenomenon would have occurred, but with different sets of indices on the u's.) 
But, in fact, analyses we have conducted indicate that it is the (even simpler) univariate 
function satisfying the relation 

SrealWl, V 2 ) = <Sreal(v) °C VVj (98) 

that fits the sample estimate of the separability function displayed in Fig. [3] extremely well. 

To substantiate this last point, in Fig. HI we show a plot of a least-squares fit of the 
normalized function shown in Fig. [3] to the function n x , which for x = | is identical to 
( |98l) . We see that the best fit does, in fact, suggest that x = ^ is the appropriate choice 
(at least, within this one-parameter family of functions). (For the least-squares fit of ([98]) 
to our sample estimate, we obtain 0.00690362, while we obtain considerably more, that 
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FIG. 4: Least-squares fit of the normalized sample estimate of the real qubit-qutrit separability 
function displayed in Fig. [3] vs. rf = {y\v%f . The minimum of the curve is in the immediate 
vicinity of x = |, at which point the measure of goodness-of-fit is 0.00690362. 

is 0.0184813, for the [inferior] fit of (j9~Tjb ) The product of the normalized function fl98|) 
with the corresponding jacobian is integrable in the real qubit-qutrit case (giving the result 
^ ~ 3.70723 ■ 10- 16 ^ 



1110124175582822400 

If we adopt the ansatz (1981) and employ the estimated value of the scaling constant for 
this function from Fig. [3j which is on the order of 3095.97, and additionally presume that the 
real qubit-qutrit HS separability probability (in line with our complex counterpart conjecture 
[immediately below] of yffg — and qubit-qubit proposals of J| and y^) is of the form ^, where 
k is some natural number, then our best estimate of this probability is ~ 0.150235, and 
of the scaling constant ^95^5^/3 ~ 3096.05. (We do not have highly extensive numerical 
estimates [only the more limited one pursued here] — as we did in the complex qubit-qutrit 
case 3] — against which to gauge this prediction, but our fits here are strongly supportive 
of these assertions. For example, our sample estimate of the separability probability can be 
expressed as 21 g q 05 -) 



B. Complex Qubit-Qutrit Case 



ly-intensive study — using 10 9 sample points — - we 



Based on our previous numerica 
have an (implicitly-given) estimate [18|, between eqs. (38) and (39)] for the complex HS 
separability probability of 0.0266891. A very well-fitting candidate for the corresponding 
exact probability is yf|g ~ 0.0266889. (The associated separable volume would, then, be 
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FIG. 5: Interpolated estimate over the unit square of the complex qubit-qutrit separability function 
Scompiexi^i-,^)-: based on 880,000 randomly generated 6x6 complex density matrices 

5698058897559008007188598937500076 ~ 2 - 05327 ' ^ ^ Aside fr ° m the StlM ^ gOodneSS-of-fit, 

we see that the numerator of the probability is equal to 32 = 2 n ~ 1 ,n = 6, while in the 
qubit-qubit case, the numerator is 8 = 2 n ~ 1 ,n = 4. Also, the denominator 1199 = 109 • 11, 
while 33 = 3 • 11. 

In line with the Dyson-indices pattern observed earlier, we investigated the possibility 
that the separability function in the complex qubit-qutrit case might be simply proportional 
to f] = v\Vi y that is, the square of its putative real counterpart, y/rj. The integral of the 
product of T) with the associated jacobian yields 5045434342262725360252343040000 ~ 1-64307- 10~ 28 . 
With our proposal above (supported by the considerable numerical evidence of |9[]) that 
the qubit-qutrit complex HS separability probability is j 3 - 2 ^, the scaling constant would be 
537472 v 7 ?- 15 ~ 1249 65 

10063956375 

Figs. [5] and [6] are the complex qubit-qutrit counterparts of the (real qubit-qutrit) Figs. [3] 
and HI Fig. [6] might be said to weakly support the proposal that the separability function 
is proportional to r\. (The numerics here are perhaps yet insufficient for our purposes. In 
addition to only so far having sampled a relatively small number of complex 6x6 density 
matrices, the sample points are now 30-dimensional in nature. For alacrity, we had simply 
used Monte Carlo methods, and not the [better-behaved/" lower-discrepancy"] quasi- Monte 
Carlo [Tezuka-Faure] methods employed in our earlier studies, in particular, in 

In 

light of the not very convincing nature of Fig. El it might be advisable to revert to the 
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FIG. 6: Least-squares fit of the normalized sample estimate of the complex qubit-qutrit separability 
function displayed in Fig. [5] vs. rf = {yiV2) x . We hypothesize that for a sufficiently large sample 
the minimum would lie at x = 1. 

Tezuka-Faure scheme, although most of the unit hypercube points generated would, then, 
be simply discarded as not meeting the criteria a density matrix must fulfill. The most 
desirable/efficient sampling scheme, it seems, if it can be effectively implemented, would be 



the one associated with correlation matrices 
points would have to be discarded. 
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36|, [37|, in which none of the generated 



XI. CONCLUDING REMARKS 



In a recent comprehensive review, it was stated that while quantum entanglement is 
"usually fragile to environment, it is robust against conceptual and mathematical tools, 
the task of which is to decipher its rich structure" 52|, abstract]. We have attempted to 
make some progress in this regard here, but considerable impediments clearly still remain 
to putting the chief conjectures of this paper on a fully rigorous basis (or disproving them 
by establishing alternative results), and in proceeding onward to higher- dimensional cases. 
(In particular, we have not yet developed a theory to predict the scaling constants — ^§§§- 



and 



207T 4 



in the full complex and real two-qubit cases, and 



537472 



and 



78848tt 



in the 



17 ^ »^ v^^^v^ v^v*,, 1006395637 5 139515^ 

full complex and real qubit-qutrit cases — for the hypothesized separability functions.) It 
would seem that applications and/or extensions of random matrix theory and, possibly, 



mathematical induction will be important — as they were in determining the total (separable 
and nonseparable) Hilbert-Schmidt volumes [lj {2], sec. 14.3] j^]. 



Let us, still, further suggest that the analytical framework and results, both theoretical 
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and numerical, presented above may lead to the development of associated formal propo- 
sitions, much in the way that the numerically-obtained (two-fold) separability-probability 



ratios reported in 



18( | led Szarek, Bengtsson and Zycskowski to establish that the set of sep- 



arable 
height 
In 



and, more generally, positive-partial-transpose) states is a convex body of constant 



20j. 



js], we have further applied the '"separability function" concept to the determination 
of the Bures (minimal monotone) metric volume of certain low-dimensional (real, complex 
and quaternionic) two-qubit states. Interestingly, we find that although the Dyson-index 
pattern is not now fully adhered too, it does come remarkably close to holding. Also, 
numerical research that we hope to shortly report, strongly indicates that in the full two- 
qubit quaternionic 2 7- dimensional Hilbert- Schmidt separable volume case, the Dyson-index 
pattern (/3 = 1,2,4) we have observed above in the two-qubit 9-dimensional real and 15- 
dimensional complex cases (Fig. [2]) is strictly maintained. 
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